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We investigate the simplest gauge theory for spontaneous R-parity breaking and its testability at 
the LHC. This theory, based on a local B-L gauge symmetry, can be considered as the simplest frame- 
work for understanding the origin of the R-parity violating interactions, giving rise to potential lepton 
number violating signals and suppressed baryon number violating operators. The full spectrum of 
the theory and the constraints coming from neutrino masses are analyzed in detail. We discuss the 
proton decay issue and the possible dark matter candidates. In order to assess the testability of the 
theory we study the properties of the new gauge boson, the neutralino decays and the main produc- 
tion channels for the charged sleptons at the LHC. We find that final states with four charged leptons, 
three of them with the same-sign, and four jets are the most striking signals for the testability of the 
lepton number violation associated with spontaneous R-parity violation at the LHC. 



I. INTRODUCTION 



The Large Hadron Collider (LHC) will hopefully soon discover the underlying theory for the TeV scale 
and might allow us to understand a more fundamental law of nature. For more than three decades the idea 
of Supersymmetry has attracted the attention of many experts in the particle physics community and the 
minimal supersymmetric extension of the Standard Model HHHl (MSSM) is still considered one of the most 
appealing candidates for the theory of particle physics at the TeV scale. It is well known that the MSSM 
provides an understanding of why the SM-like Higgs boson is light, contains a cold dark matter candidate, 
allows for the unification of the gauge couplings and allows for the mechanism of electroweak baryogenesis 
to explain the baryon asymmetry in the universe. 

There are several open issues in the MSSM, one of them being the origin of the discrete symmetry R- 
parity |@J|5l. This symmetry plays a major role in the MSSM and it is defined as R = (— l) 3 (- B_i )+ 2S ', 
where B, L and S stand for baryon number, lepton number and spin, respectively. In many MSSM studies it 
is assumed that R-parity is conserved or explicitly broken without understanding the origin of this symmetry. 
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However, the fate of R-parity is crucial for the discovery of supersymmetry since, as is well known, R-parity 
conservation give rise to channels with multi-jets, multi-leptons and missing energy at the LHC, while 
signatures of broken R-parity are multi-leptons, multi-jets, and missing energy due to the SM neutrinos 
only. 

The simplest and most elegant framework for the origin of R-parity is based on local B-L symmetry. 
This connection was explored for the first time in Ref. and in Ref. [7] a simpler scenario was studied. 
See also Ref. [ 8 ] for a complete discussion of how to gauge R-parity. 1 Recently, we have investigated the 
simplest B-L models in Refs. |[T4] - [T7 l and found the following main result 

The simplest theories based on local B-L make the following prediction: 
R-parity must be spontaneously broken at the TeV scale and 
one expects to observe lepton number violation at the LHC ! 

In this letter we study in detail the theory proposed in Ref. lTT5l which can be considered as the simplest 
gauge theory for R-parity violation. In this context the only way to break local B-L and obtain the MSSM 
after symmetry breaking is to give a vacuum expectation value to one of the right-handed sneutrino required 
by anomaly cancellation. One of the most important features of this theory is that the B-L and R-parity 
breaking scales are determined by the soft supersymmetric breaking scale. This idea was studied for the 
first time in Ref. lfl4l which defined the simplest left-right symmetric model. 

We review the theory and symmetry breaking mechanism in Sections II and III. The full spectrum of 
the theory |[T5l is discussed in Section IV and the constraints coming from neutrino masses in Section VI. 
We discuss the proton decay issue and the possible dark matter candidates in Section V. In order to under- 
stand the testability of the theory we study the properties of the new gauge boson, the neutralino decays in 
Section VII and in Section VIII the main production channels for the charged sleptons at the Large Hadron 
Collider. We find that the channels with four charged leptons (three with the same electric charge) and four 
jets give us the most striking signals for the testability of lepton number violation at the LHC. 

II. THE MINIMAL GAUGE THEORY FOR SPONTANEOUS R-PARITY VIOLATION 

The simplest gauge theory for spontaneous R-parity breaking was proposed in Ref. |[T5l . In this context 
one can understand dynamically the origin of the R-parity violating terms in the MSSM. Here we discuss 
the structure of the theory and the full spectrum. 



' It is important to mention that the breaking of B-L in the context of the MSSM was studied for the first time in Ref. 1 9 1 . See also 
Refs. 1 10-13 1 for the study of R-parity in other models. 



Gauge Symmetry and Matter Fields: This theory is based on the gauge group 

SU(S) C ® SU(2) L U(1) Y ® U(1) B -l, 
and the different matter chiral superfields are given by 

Q= (^j ~ (2,1/3,1/3), L= ~ (2,-1,-1), (1) 

n c ~ (1,-4/3,-1/3), cf ~ (1,2/3,-1/3), e c ~ (1,2,1). (2) 

In order to cancel the B — L anomalies one introduces three chiral superfields for the right-handed 
neutrinos: 

v c ~ (1,0,1). (3) 
Higgs Sector. The Higgs sector is composed of two Higgs chiral superfields as in the MSSM 

H u =\ H ") ~ (2,1,0), H d = I Hd I ~ (2,-1,0). (4) 

Interactions : With this field content the superpotential reads as 

W BL = Wmssm + Y u L T ia 2 H u C , (5) 

where 

Wmssm = Y u Q T ia 2 H u u c + Y d Q T ia 2 H d d c + Y e L T ia 2 H d e c + fi ia 2 H d . 



(6) 



In addition to the superpotential, the model is also specified by the soft terms: 



t t 2 i~c|2 . 2 

V so ft = m^v | + m~ L 



L 



2 + ml |e c | 2 + m 2 Hu \H u \ 2 + m 2 Hd \H d \ 2 + (±M BL B'B' 



+ A u L T ia 2 H u D c + Bfi B T U ia 2 H d + h.c.) + V s %f M , (7) 

where the terms not shown here correspond to terms in the soft MSSM potential. 

Since we have a new gauge symmetry in the theory we need to modify the kinetic terms for all 
MSSM matter superfields, and include the kinetic term for right-handed neutrino superfields 

C K in{y c ) = j d 2 6d 2 6 {u c )^e 9BLVBL u c . (8) 

Here Vbl is the B-L vector superfield. Using these interactions we can study the full spectrum of the 
theory. 
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III. ELECTROWEAK AND B-L SYMMETRY BREAKING 

As in the MSSM, electroweak symmetry is broken by the vevs of H® and H®, while U (1) b-l is broken 
due to the vev of right4ianded sneutrinos. Notice that this is the only field which can break local B — L and 
give mass to the new neutral gauge boson in the theory. Therefore, the theory predicts spontaneous R-parity 
violation. It is important to mention that the B — L and R-parity breaking scales are determined by the soft 
supersymmetric breaking scale, and one must expect lepton number violation at the LHC. 

The neutral fields are defined as 



-j= {vd + h d ) + ~^ A d, 



and the relevant scalar potential reads as 



V 
V F 

Vd 

Vsoft 



V F + V D + Vsoft, 

H 2 |^| 2 + i-MflS + W^I 2 + EttJ^lW + ElW'lW 

i 3 



l^l 2 -I^T-EN 2 ) + 9 -f 



Bi*fi 2 -ifci 2 



(5f) f 



l / \ z / 



(9) 
(10) 

(11) 
(12) 



(13) 
(14) 



*i| 2 )J , (15) 
r0 u - BfiH^H® + h.c.) . 



(16) 



Using the above scalar potential and assuming that all parameters are real we can find the minimization 
conditions 



/, 2 + -Y?tf R Y?v k R + 2 4^'^n fej + 
+ ^i< j v R - Bfiv d = 0, 



Vd 



/i 



2 (01 + S2) / 2 



pv i L Yj*v' R - Bixv u = 0, 



(17) 
(18) 
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2 \ V L Y v V R V L Y v 

' [(mlc) ki + (mlc) ik ] + —victfvv = 



9sl 



J R 



V2' 



(19) 



<y^4^ m < 



n - 44) 4 



(44 - 44) 4 + J 4 + 



+ ^4^>« 



o. 



(20) 



In order to have phenomenological allowed solutions the have to be small, and the v l R have to be much 
larger than v u ,Vd and v\. Up to negligibly small terms 2 the right-handed sneutrino acquire a vev in only 
one family. A possible solution and the one used throughout this paper is v R = (0, 0, vr). In this case: 

(21) 
(22) 



4 - 


. 8(m? c ) 33 






i > 
9bl 




4 - 






{\xv d Y^ - afv u ) 


' V2 


{m 2 L )kk 


(9l+s|)c,,2 ,,2^| 9Il ? ,2" 
8 °d) 8 U i? 



Notice that the minimization conditions for v u , Eq. ( (T7| ), and tj^, Eq. ( |18) are not greatly altered from their 
MSSM equivalents since the extra terms are very small. 



A. Radiative Symmetry Breaking 

In the MSSM, the large top Yukawa coupling drives the up-type soft Higgs mass squared parameter to 
negative values for generic boundary conditions leading to radiative electroweak symmetry breaking [18]; 
a celebrated success of the MSSM. A valid question is then if the same success is possible in achieving a 
tachyonic right-handed sneutrino mass in this B — L model as required by Eq. ( |2T[ ). Unfortunately, this is 
not possible through a large Yukawa coupling since the Yukawa couplings of the right-handed neutrino are 
all dictated to be small by neutrino masses. However, there is an alternate possibility whereby a positive 
mass squared parameter for the right-handed sneutrino at the high scale will run to a tachyonic value at the 
low scale. This is due to the presence of the so-called S'-term (due to D-term contributions to the RGE) in 
the soft mass RGE, as discussed for this B — L model in (T9J-ED- A short outline of the mechanism follows. 

The RGE for the right-handed sneutrino soft mass squared parameter is 

16vr 2 ^ = -3gl L \M BL f + -/ BL S BLl (23) 

The size would go as y »^« v l < 1(T 10 . The maximum values for Y u and Vl are about 10 6 and 10 2 GeV, respectively, see 
Fig.|4] 
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with 

Sbl = Tr (2mJ - m? c - m| c - 2m| + mfc + m~ c ) , (24) 

where the trace is over the three generations of the fermions and the soft mass parameters in the trace are 
for the squark doublet, right-handed up squark, right-handed down squark, slepton doublet, right-handed 
charged slepton and right-handed sneutrino, respectively. The gaugino mass term always drives the sneu- 
trino mass parameter positive at the low scale but if the overall sign of the S-term is positive, it could lead to 
the opposite effect. Such an effect would require a non-zero S'-term at the high scale, which is not possible 
if the soft masses are universal across the generations of each flavor. An example of a suitable boundary 
condition with minimal variation from the popular MSUGRA Ansatz is universal boundary conditions for 
all sfermions except for the right-handed sneutrinos, which might have the boundary conditions 

m lf = m \ = p ml, m|c = Qml, (25) 
where mo is the universal mass and P > 1 and Q < 1. The boundary condition for Sbl is then 

S BL = (2P + Q-l)m 2 , (26) 

having the necessary sign to contribute negatively to the sneutrino soft mass parameter as it is evolved from 
the high scale down. The necessary sizes of P and Q depend on the size of the gaugino mass parameter 
which has the opposite effect, see IT9H2TJ for more details. So while, the traditional radiative symmetry 
breaking from universal boundary conditions is not possible in these models, it is possible to radiatively 
break B — L through this S-term starting from a positive value. For the implementation of the radiative 
mechanism in the non-minimal model see Ref. Il22ll . 



IV. MASS SPECTRUM AND LEPTON NUMBER VIOLATION 
A. R-Parity Violating Interactions 

After symmetry breaking lepton number is spontaneously broken in the form of bilinear R-parity vio- 
lating interactions. There are no trilinear R-parity violating interactions at the renormalizable level. These 
bilinear interactions mix the leptons with the Higgsinos and gauginos: 

\gBLV R (v c 3 B), \g2vi(yiW°), -±=g 2 vi(e t W + ), 
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The first term is new and is the only term not suppressed by neutrino masses. The fifth term corresponds to 
the so-called e term, and second, third and fourth terms are small but important for the decay of neutralinos 
and charginos. See Section V for the discussion of the baryon number violating operators. 

There are also lepton number violating interactions coming from the soft terms and the B-L D-term. 
From V so ft one gets 

Af^Lf ia 2 H u , 

while from the D-term one finds 

9 2 blVrv c ($\<i- ■ 
As one can expect these terms are important to understand the scalar sector of the theory. 



B. Mass Spectrum 

Gauge Boson: The neutral gauge boson associated to the B — L gauge group is Zbl- Using the covariant 
derivate for the right-handed sneutrinos, = d^v c + i^^B^D , the mass term for Zbl is: 

Mz BL = 9 -^vr- (27) 
Now, using the experimental collider constraint ||23ll : 

> 3 TeV, (28) 



Mz BL 



9BL 

and Eq. (21) one finds the condition 

|(m 5 c )33 | > 2.12 g BL TeV. (29) 

Then, if gsL = 0.1 the soft mass above has to be larger than 200 GeV. This condition can be easily satisfied 
without assuming a very heavy spectrum for the supersymmetric particles. 

Neutralinos and Neutrinos: As in any supersymmetric theory where 72-parity is broken all the fermions 
with the same quantum numbers mix and form physical states which are linear combinations of the original 
fields. The neutralinos in this theory are a linear combination of the fields, ( i/j, v?, B' , B, W°, H®, H®). 
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Then, the neutralino mass matrix is given by 



M 



N 






72 Yu Vu 


1 ? 

~29BL Vl 


1 4 

~29i vl 


1 i 

292 Vl 





V2 Yv v r 


1 v l i D 

V2 V Vu 





hBL VP R 











V2 Yu V L 


1 ■>' 
~29BL V* L 


1 7 
2.9BL V J R 


M BL 














-hi vi 








Mi 





-hwd 


55i^« 


h? V L 











M 2 


h^ v d 


-\giv u 













l92 v d 







J_ yij 3 
\V2 V R 


r " U ^ 





551 


-\92 v u 








(30) 



/ 

We have discussed above that only one right-handed sneutrinos get a vev, v l R = (0, 0, vr). Now, integrating 
out the neutralinos one can find the mass matrix for the light neutrinos. In this case one has three active 
neutrinos and two sterile neutrinos, and the mass matrix in the basis (i/ e , v^, u T , z/g, u^) is given by 



where 



My 



A 



V 

m 3 



A vlv L + B 



B 



Y?vl 



\V2vr 



+ 



V2/1V, 



dVR 



m J 



+ CY?Yj A ^v u Yy 



C 



o 



2x2 



( 2M BL v 2 u 
V 9bl v r 



+ 



v d v R 
TO 3 



(3D 



(32) 



m 3 



4 [nv u v d {g\M 2 + glM x ) - 2M 1 M 2 ^ 2 



(33) 



g\M2 + g\M x 

Here a and /3 take only the values 1 and 2. From the experimental limits on active neutrino masses we 
obtain {Yv)ia ^ 10~ 12 . This can be compared to {Y v )^ < 1CT 5 , which is less constrained because of the 
TeV scale seesaw suppression. It has been pointed out in Ref. [24] [25j (and earlier in a different context Q) 
that this theory predicts the existence of two light sterile neutrinos which are degenerate or lighter than the 
active neutrinos, a so-called 3 + 2 neutrino model. A sample possible spectrum is displayed in Fig. [I] 

Recently, it has been shown in ||26l that precision cosmology and big-bang nucleosynthesis mildly favor 
extra radiation in the universe beyond photons and ordinary neutrinos, lending support to the existence of 
sub-eV sterile neutrinos. 

Charginos and Charged Leptons: In this theory the chargino mass matrix, in the basis ^e^, W + , H^J 
and (ej, W~, H7 ), is given by 



M c 
Ml 



(34) 
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Possible spectra: 
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^2 
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FIG. 1: Sample spectra for neutrino masses in the normal and inverted hierarchies. 



with 



' V2 e d 



M c 





M 2 



1 Y e ij ^ j -\ 
i 



V2 Ye V L 



T2 92V d 



Squarks and Sleptons: In the sfermion sector, the mass matrices .M ~, and Ai - for squarks, 
charged sleptons, in the basis ( /, f c * ) , are given by 



(35) 



and Ml for 



Ml 



3 5 s^) M§ c 2/ g + \D BL (a„ v u -Y u fi v d ) 

^ (a« u u - y u fi v d ) m\ c + ml + §Mf c 2/3 s^, - 



(36) 



Ml 



m l+ m d~ {\ ~ \ *w) M Z C2/3 + \D BL 



7= (Irf fiv u - a d v d ) 



v/2 



^ (*d fJ-v u - a d v d ) 



'I ~ s w0 M i c 2/3 - -^bl -75 (Y e ixv u -a e v d ) 



7= (y e fj, v u - a e v d ) 



v/2 



where c 2 /j = cos 2/3, s\y = sin and 



Del 



1 



2 2 



1 



-MS 



(37) 



(38) 



(39) 



m u , m d and m e are the respective fermion masses and a u , a d and a e are the trilinear a-terms corresponding 
to the Yukawa couplings Y u , Y d and Y e . Typically, it is assumed that substantial left-right mixing occurs 
only in the third generation. Regardless, the physical states are related to the gauge states by 

fl\ _ I cos6 / sin9 A I f 
J 2 [-smOf cose? \r 



(40) 
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The masses in the sneutrino sector are given by 



Ml 
Ml 
Ml 



ml 



Li + - Ml cos 2/3 -~Mt BLl 



M 



m 



2 i 



D BL , 



(41) 
(42) 
(43) 



and a = 1..2. For simplicity we listed the mass matrices in the limit v l L , a v , Y v — > 0. For the most general 
expressions see Appendix A. It is important to mention that all sfermion masses are modified due to the 
existence of the B-L D-term. 

In order to understand the properties of the spectrum we assume a simplified spectrum for the super- 
partners. In the case of the sfermions we will assume for simplicity the same value for all soft masses. In 
this case if we neglect the left-right mixing the full spectrum of sfermions will be defined by Mstjsy (the 
universal soft supersymmetry breaking mass), tan /3 and the mass of the B-L gauge boson. Using this sim- 
plified spectrum we show in Fig.l the values for the sfermion masses for different values of Mz BL - Notice 
that the condition that left-handed slepton masses have to be positive impose a bound on the Mz BL for a 
given value of Msusy, i-e. Mz BL < \f2Mp In this way, we can appreciate that the spectrum can be very 
constrained. As it is well-known, in the general case one cannot predict the spectrum since the soft masses 
are unknown. 
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FIG. 2: Spectrum for sfermion masses assuming the same value for all soft masses, Msusy — 1 TeV and tan/3 = 5. 
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V. NUCLEON STABILITY AND DARK MATTER 

It is well-known that when R-parity is broken in a supersymmetric theory one has to understand the 
possible constraints coming from proton decay E71 . In the MSSM one has several interactions which could 
mediate proton decay at tree level and one-loop level. At the renormalizable level one has the lepton and 
baryon number violating interactions 

W RpV = eLH u + XLLe c + XQLd c + \"u c d c d c , (44) 

which are not allowed in our theory before B-L breaking, and in general one has the dimension five operators 

W RpC = j-LLH u H u + ^QQQL + ^u c d c u c e c + ^u c d c d c u c . (45) 

Notice that the first term in the above equation is not allowed in our theory, but the last terms can mediate 
proton decay. The operators QQQL and u c d c u c e c mediate proton decay at one-loop level and typically 
the scale A should be larger than 10 17 GeV in order to satisfy the experimental bounds on proton decay. 
For a detailed discussion see Ref. E71 . Once B-L is broken by the vev of the right-handed sneutrinos one 
finds new contributions to proton decay at tree level. From the Yukawa coupling Y u LH u v c one gets the 
lepton number violating interaction LH U and from the last term, u c d c d c u c , in the above equation one gets 
the interaction u°d c s c . Using these interactions and integrating out the neutralinos and squarks we find the 
following constraint 

from the channel p — > K + v. Then, assuming that A^ 23 ~ 1, and Y"J 3 ~ 1CT 6 (see Fig. [4]) one gets 
A > 10 17 GeV. This constrain is similar to the one we get from the dimension five operators. Therefore, 
one can say that if the above couplings are order one the cutoff of the theory has to be large. Also one can 
think about possible scenarios where the couplings are small, see for example ll28l . 

At first glance, finding a dark matter in .R-parity violating theories seems hopeless. But while the 
traditional neutralino LSP case is no longer valid, the situation is not lost. As first discussed in 112911 , such 
models can have an unstable LSP gravitino, with a lifetime longer than the age of the universe. The strong 
suppression on its lifetime is due to both the planck mass (Mp) suppression associated with its interaction 
strength and bilinear i?-parity violation which is small due to neutrino masses and must facilitate the decay 
of the LSP. In the mass insertion approximation, this can be understood as the gravitino going into a photon 
and neutralino which then has some small mixing with the neutrino due to i?-parity violation (m x v ), thereby 
allowing G — > as in Fig. 3. Adopting approximations made in [29], the lifetime for the gravitino 
decaying into a photon and neutrino (in years) is about 
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FIG. 3: Gravitino decay into a photon and a neutrino. 



which for appropriate values of the gravitino mass leads to a long enough life time. Unlike in imparity 
conserving models with a gravitino LSP, there are no issues with big bang nucleosynthesis from slow NLSP 
decay since the NLSP decays more promptly through i?-parity violating interactions. Several interesting 
studies have been conducted on the signatures and constraints of unstable gravitino dark matter, see for 
example 



VI. EXPERIMENTAL CONSTRAINTS 



In order to understand the different lepton number violating decays in the theory we need to understand 
which are the main constraints from neutrino experiments. Today, we know well the numerical values for 
two of the neutrino mixings and the mass squared differences. The neutrino mixing matrix Vpmns is 
defined as 



/ 



VpMNS 



\ 



(48) 



Cl2 C13 Ci 3 Sl 2 Sl3 

-C23 S12 - Ci 2 Si3 S 2 3 C\ 2 C 2 3 ~ «12 «13 «23 C13 S 2 3 
\ S12 S23 - C12 C23 S13 -C12 S23 ~ C23 «12 «13 C13 C23 J 

where c-ij = cos 6*^ and Sij = sin 0y with < 0y < n/2. The physical neutrino masses are contained in 
m v = diag(m Ul , m U2 , m V3 ). As it is well-known, there are two possible neutrino spectra: 



Normal Hierarchy (NH): m Ul , m V2 
Inverted Hierarchy (IH): m Ul - 



m'f,. + Am 2 , , in,,. 



m 2 3 + |Am 2 



311' 



2D '"va 
m U2 = 



311' 



m 2 ^ + | Am 2 



(49) 



where OT 



7.27 x 10~ 5 eV 2 < Am^ < 8.03 x 1(T 5 eV 2 



2.17 x 1CT 3 eV 2 < |Am§i| < 2.54 x 1CT 3 eV 2 , 



(50) 
(51) 
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are the solar and atmospheric mass squared differences, respectively. 

In order to understand the allowed values for the vevs of the left-handed sneutrinos and the Dirac Yukawa 
couplings we assume for simplicity that the off-diagonal block matrices in Eq. ( |3~Tj ) are zero, hence decou- 
pling the two light sterile neutrinos from the active ones. In this case, all neutrino masses and mixing 



originate from the upper-left block matrix in Eq. (31), which we label m u . The flavor pattern, and hence the 
rank of this matrix, dictates that one neutrino will be massless. This matrix is diagonalized by the PMNS 
matrix 

m u = Vpmns M v Vpmns, (52) 



where m v = diag(0, \J Am^, y|Am|j) in the Normal Hierarchy and m„ 
diag(Y / |Am3 1 |, y^Am^l + Am^, 0) in the Inverted Hierarchy. This yields a system of six equa- 
tions quadratic in the vevs of the right-handed sneutrinos and Yukawa couplings, although solving for these 
is not the most efficient way to proceed. Instead, notice that the product above yields the following six 
terms 



(53) 



yj = yi3 yij 
1 — 1 v ''PMNS- 

Normal Hierarchy 

In the Normal Hierarchy one obtains the following six equations 

AV? + 2BV X Y X + CY? = 0, (54) 

AViV 2 + B (ViY 2 + V 2 Y l ) + CY t Y 2 = 0, (55) 

AViV 3 + B (ViY 3 + V 3 Yi) + CY X Y 3 = 0, (56) 

AV 2 V 3 + B (V 2 Y 3 + V 3 Y 2 ) + CY 2 Y 3 = 0, (57) 

AV 2 2 + 2BV 2 Y 2 + CY 2 2 = m 2 , (58) 

AVi + 2BV 3 Y 3 + CY 2 = m 3 . (59) 

In order for these equations to be consistent, V% = Y\ = 0. While this condition represents some 
fine-tuning between parameters, it is a result of the simplifying assumption that the sterile and active 
light states decouple. In a more general scenario, this condition would not be necessary. 

The remaining system of equations, the last three, is underdetermined with three equations and four 
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unknowns. Solving with respect to Y 3 yields 

Y 2 = ei 

V 3 = - 
V 2 = - 



-Y 3 2 R - Am 3 



7Tt3 V R 

BY 3 + e 3y /Y*R + Am 3 
A 

BY 2 + e 2 ^Y 2 R + Am 2 



with 



R = B 2 — AC, ei = ±1, e 2 = ±1, e 3 



ei RY 3 



e 2 7i^72' 



(60) 

(61) 
(62) 

(63) 



Inverting Eq. (53) one translates these solutions to the variables of interest. The result being that 
specifying the SUSY spectrum and B-L parameters {Mz BL and qbl) as well as Y 3 specifies all the 
values for v l L and Y"J 3 . 



Inverted Hierarchy 

In the case of the inverted spectrum one can use the same procedure. However, one needs to make the 
following replacements: m\ <-> m 3 , Y\ Y 3 , V\ V 3 . In this way when we solve the equations 
for Vi and Y{, one obtains V 3 = 0,Y 3 = and the solutions above where we have made the previous 
substitutions. 



Parameter 


Range 


Mi 


100 - 1200 GeV 


M 2 


100 - 1200 GeV 


ImI 


100 - 1200 GeV 


tan j3 


3-50 


\Y 3 \ 


10~ 7 - 10~ 5 


M BL 


100 - 1200 GeV 


M Zbl 


1000 GeV 


9bl 


0.33 



TABLE I: Parameters and ranges scanned for the plots in Figs. HI p- 13 



In Fig. |4] we show the allowed values for the vevs of the left-handed sneutrinos and the Dirac Yukawa 
couplings for a scan over the parameters listed in Table [i] As one can appreciate the allowed values for v l L 
are in the range (10~ 2 — 10) MeV while the Yukawa couplings change between 10~ 7 and 10~ 5 . Now, using 
these results we are ready to discuss all R-parity violating decays. 
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Normal Hierarchy 
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Inverted Hierarchy 
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(e) 




Kl (MeV) 
(f) 



FIG. 4: Allowed values for the v l L versus (Y„)i3 in agreement with the neutrino masses and mixings constraints in 
the NH (IH) in a,c and e (b,d and f) for a scan over the parameters listed in Table [I] 
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VII. LEPTON NUMBER VIOLATION AND DECAYS 

At this point, the relevant pieces of this model have be laid out and the question of interesting signals 
can now be tackled. Since lepton number is violated, same-sign dileptons and multijet are possible final 
states. Such signatures are interesting since they have no SM background. However, the final states depend 
critically on the nature of the LSP and since R-parity is violated the possibilities are more numerous than 
normal, i.e. colored and charged fields. 

These possibilities are briefly discuss in Appendix[Dj We find that the most clear single of lepton number 
violation (and therefore the most interesting for us) results from the decays of a neutralino LSP through the 
process (see Fig [5]) 3 : 

pp -> 7*, Z*, Z* BL -> e*e, -> ef ef ej ef 4j. 

In order to quantify this signal we will continue by investigating the decays of the Zbl gauge boson, 
the charged slepton and the neutralino and in the next section, the production mechanism for the charged 
sleptons at the LHC. 




FIG. 5: Topology of the signals with multi-leptons. 



3 Throughout this paper, shorthand such as pp — > e*e represents the process pp — > e* e + X, where the activity associated with X 
has low transverse momentum and is not associated with the relevant physics of interest. Alternatively, this notation represents 
all possible production methods of e*e from the partons inside the proton taking their respective parton distribution functions 
into account. 
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A. B-L Gauge Boson Decays 



The Zbl boson can decay into a pair of charged fermions, light neutrinos, and into two sfermions. The 
partial widths for the decay into particles Pi, P2 of masses mi, rrt2 are given by 



T(Z BL ^P 1 P 2) 



16ttM Zp 



\M(Z BL P!P 2 )| 



(mi + m 2 ) 2 

— M 

^BL 



(mi — m<iy 



M 



7 2 



(64) 



where the squared matrix elements for specific final states are in Appendix [B] The Zbl branching ratios 
are plotted in Fig. [6] for a fixed soft universal mass for all sfermions, Msjjsy = 1 TeV, versus the Zbl 
mass. Decay channels into SUSY particles only open up for a Zbl mass of around 1.2 TeV and at much 
1.4 TeV, the sleptons become tachyonic. As can be appreciated from Fig. [2], tachyonic sleptons are reached 
before decay channels into the squarks can open. The branching ratios for the sleptons are divided up 
into sneutrinos, smuons plus selectrons and staus in anticipation of the associated signals. However, each 
individual slepton pair has the same branching ratio, about 2.5% at mz BL = 1.4 TeV. In Fig. 6 we show 
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FIG. 6: Zbl branching ratios versus Zbl mass for a universal soft mass for all sfermions, Msusy = 1 TeV. Here 
the subscript one refers to the lightest eigenstate in each family and this case corresponds to the purely left-handed 
slepton (zero mixing angle is assumed). Given the universal soft mass for the sfermions, only the left-handed slepton 
channels can be open. Right-handed squark channels can open for larger values of Mz BL but only at the cost of 
unphysical tachyonic slepton masses. 



the predictions for the total decay width and the invisible decays of the Zbl gauge boson. Since the new 
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gauge boson can decay into five light neutrinos the invisible decay can be large, a few GeV when the mass 
is above 1 TeV. These properties of the Zbl are very important in order to discover this theory at the LHC. 
In summary, one can say that the new neutral gauge boson is B-L like with branching ratios 

Br(Z BL -> e+er) ~ 40%, Br(Z BL -> vv) ~ 35%, Br(Z BL -> jj) ~ 20%, and Br(Z BL -> tt) ~ 5%, 

since the branching ratios for the SUSY decays are small and the invisible decay can be large. For example 
when Mz BL = 1.2 TeV the invisible decay width is Tz BL (invisible) ~ 3 GeV. 




200 400 600 800 1000 1200 1400 

M ZBL (GeV) 

FIG. 7: The total decay width, T(Zbl), versus Zbl mass for a universal soft mass for all sfermions, Msusy = 1 
TeV. 



B. Charged Slepton Decays 

The leading decay channels for the charged left-handed sleptons, ef , are the decays into neutralinos and 
charginos 

ef^ e tx°a, ef^vx% ef^fyW*, 

where i,j are lepton generational indices, a labels the neutralinos from lightest to heaviest and A labels the 
charginos from lightest to heaviest. In addition to these decay modes there are various imparity violating 
decays which only dominate when the slepton is the LSP. The last channel above usually involves an off- 
shell product particle (a three-body decay) and is therefore suppressed. The decay widths for the remaining 
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two relevant channels are given by 




9BLN 2a + 9i N 3a + g 2 N4a\ + 



2V2m e 

Vd 



5(i 




m|o' 
ml 



v iXa) 



16vr 



rat 



1AI 



1 



(65) 



(66) 



where iV a fr diagonalizes the neutralino mass matrix and in the chargino sector one has V~XV + = 
diag (m^ 1; m^ 2 ) . Due to the presence of the right-handed neutrino and B', position 2, 3, 4 and 5 in the 
N refer to B', B, W and, H^, respectively. While it is hard to make predictions for the branching ratios of 
the charged sleptons decaying into a charged lepton and LSP without knowing the details of the spectrum, 
we briefly outline the best case scenario. Our final results will be given both with this best case scenario in 
mind and with arbitrary Br(e« — > ejXi)- 

For a mostly bino LSP, it is possible that charged selectrons decay one hundred percent into the LSP 
since the charginos and other neutralinos could be heavier. For a wino LSP however, the lightest chargino 



channel is very likely to be open. Because of the factor of 2 difference between Eq. (65 1, and Eq.(66l if all 
other neutralinos and charginos are kinematically disallowed, one can expect 

T(ef -+ efxl) 1 



(67) 



meaning a 33% branching ratio for a charged slepton into a wino LSP. Meanwhile, the left-handed selectron 
does not couple to the charged Higgsino. Therefore, in the same limit where all other neutralinos and 
charginos are out of kinematic range, the charged sleptons decays one hundred percent into the LSP for a 
Higgsino LSP. 



C. Neutralino Decays 

The leading decay channels for the lightest neutralino, x?, include 

xl^e±W*, Xi^ViZ, xi^Hhk, Xi^efH*. (68) 

The amplitude for the two first channels are proportional to the mixing between the leptons and neutralinos, 
while the last one is proportional to the Dirac-like Yukawa terms. While decays to all the MSSM Higgses 
are possible, typically, only the lightest MSSM Higgs, h (k = 1), is light enough for the scenario we 
consider here and so we will only take it into account. 
A naive estimation of the decay width yields 

r (x?)~ |HK x | 2 M x , (69) 
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where V xv is the mixing between the neutralino and neutrino which is proportional to y / m u /M x . Assuming 
that m v < 0.1 eV the decay length one finds L(x°i) > 0.6 mm. Therefore, even without making a 
detailed analysis of the decays of the lightest neutralino one expects signals with lepton number violation 
and displaced vertices in part of the parameter space. For a recent analysis of the neutralino decays in 
R-parity violating models see Ref. 



The specific decay width expressions are 



= T(x a efWf) = ^W^ml ( 1 - ^ ) , (70) 
T e,w T = ma ^ efw ^ = A-\ Via fm Xa (l-^f) , (71) 
1^ - r(Xa -+ ViZ L ) = ^W^ml [l-^f), (72) 





r^Erfc^^ T ) = ^|v,|^,,|i-^| . (73) 

T^ ee T(x a -+ nh) = ^r\V ia \ 2 cos 2 a m| o ( 1 - ££- I . (74) 

Here a is the mixing angle in the Higgs sector and in the decoupling limit, M\ » M|, which we assume: 
cos a = sin /?. The index i indicates the generation of lepton and a the neutralino with a = 6 is the 
heaviest and a = 1 is the lightest. These expressions depend on the mixing between the light neutrinos and 
the neutralinos, Vi a , which is derived in Appendix [C] Of course, only the decays of the LSP are relevant 
since the decays of the other neutralinos will be dominated by 7?-parity conserving decays, hence a = 1 
for our purposes. Table [I] displays the values of interest for a specific point in parameter space to gain an 
appreciation for possible values. Notice that in these scenarios the branching ratios for the channels with 
charged leptons can be large. 

In Figs. [8 10 are the decay lengths versus LSP mass resulting from a scan over all the possible values of 
ei and and over the parameters and ranges specified in Table [I] The points are divided according to the 
largest component of the LSP and the neutrino hierarchy with a dominantly bino, wino and Higgsino LSP 
in the NH shown in (a) and for an IH in (b), respectively. The relevant decay lengths can be understood 
by studying the mixings in Eq. ( |30| ). Since the higgsino-neutrino decay strength is the largest, ~ Y u vr, 
the Higgsino LSP has the shortest decay length. It is followed by the wino LSP with mixing ~ Q2Vl and 
finally the bino with coupling ~ giV£ and therefore the largest possible decay lengths. Displaced vertices 
associated with the lifetime of the LSP will only be discernible in a very limited part of the parameter space. 

The LSP branching ratios into the various possible channels versus the LSP mass are displayed in 
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Hierarchy 


LSP Personality 


Bino 


Wino 


Higgsino 




Decay Length 


1.1 mm 


0.03 mm 


1 x 10 -4 mm 




2Br(x? -> e~W+) 


4% 


2% 


13% 




2Br( X ? -> li- W+) 


26% 


12% 


27% 


NH 


2Br(x? ^t-W+) 


61 % 


54% 


30% 




Br( X ? -> 


10% 


29% 


28% 




Br( X ? -> i>fc) 


0% 


3% 


1 % 




Decay Length 


0.6 mm 


0.01 mm 


1 x 10~ 5 mm 




2Br(x? -^erW+) 


17% 


3% 


25% 




2Br(x? 


36% 


32% 


19% 


IH 


2Br(x? -> t-W+) 


38% 


34% 


26% 




Br(x? -> 


10% 


29% 


28% 




Br(x? -> vh) 


0% 


3% 


1 % 



TABLE II: Values of interest for a sample point in parameter space: ei = £2 = 1,1a = Y\ = 10 6 ,Alz BL = Mg, = 
1 TeV, tan/3 = 5 and m h = 125 GeV. Here M X ,M 2 and fi are 100 GeV, 500 GeV and 500 GeV for the bino LSP 
case, 500 GeV, 150 GeV and 500 GeV for the wino LSP case and 500 GeV, 500 GeV and 100 GeV for the Higgsino 
LSP case respectively. 
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FIG. 8: Decay length in millimeters versus LSP mass for a dominantly bino LSP in (a) for a NH and in (b) for an IH. 
Parameters are scanned according to the ranges specified in Tableland over all values of ej and e 2 . 



700 800 



Figs. TT][T3 scanning over the parameters in Table [I] and plotting the dominantly bino, wino and Higgsino 
LSP in (a) for a NH and in (b) for an IH. The lack of variance with scanned parameters displayed in the 



22 





IH: Wino LSP 
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FIG. 9: Decay length in millimeters versus LSP mass for a dominantly wino LSP in (a) for a NH and in (b) for an 
IH. Parameters are scanned according to the ranges specified in Tableland over all values of e\ and e^- 



IH: Higgsino LSP 




100 200 300 400 500 600 700 
LSP Mass (GeV) 
(a) 



800 



100 200 300 400 500 600 700 800 
LSP Mass (GeV) 
(b) 



FIG. 10: Decay length in millimeters versus LSP mass for a dominantly Higgsino LSP in (a) for a NH and in (b) for 
an IH. Parameters are scanned according to the ranges specified in Tableland over all values of t\ and 62. 

v Z and v h channels are due to the sum over all three flavors of neutrinos and also exists for the sum over 
the three charged lepton plus channels which total about 50% (or more if the other channels hadn't 
fully turned on yet). Although it is not obvious from the Figs. TT|[T3 the branching ratio to the electron 
channel is always smaller then either the pL^W^ or the t^W^ in the NH. Since we now know the 
properties of the neutralinos and selectrons decays we are ready to study the production channels. 
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FIG. 1 1 : LSP branching ratios versus LSP mass for a dominantly bino LSP in (a) for a NH and in (b) for an IH. 
Parameters are scanned according to the ranges specified in Tableland over all values of ej and e 2 . 




FIG. 12: LSP branching ratios versus LSP mass for a dominantly wino LSP in (a) for a NH and in (b) for an IH. 
Parameters are scanned according to the ranges specified in Tableland over all values of t\ and t2- 

VIII. PRODUCTION MECHANISMS AND SIGNALS 

The lepton number violating signal discussed in the beginning of the previous section proceeds from the 
pair production of charged sleptons. While the MSSM contributions to this production drop rapidly with 
charged slepton mass (since that mass must be above the Z threshold), in this model further contributions 
due to the Z' resonance can significantly increase the cross section. These enhancement are discussed first 
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FIG. 13: LSP branching ratios versus LSP mass for a dominantly Higgsino LSP in (a) for a NH and in (b) for an IH. 
Parameters are scanned according to the ranges specified in Tableland over all values of ej and e 2 . 

followed by a study of the expected number of events at the LHC. 



A. Sleptons Production Mechanisms 

The main production channel for the charged sleptons is through the photon, the Z gauge boson and the 
Zbl boson 

q(Pi)q(P2) -> J,Z*,Z* BL -)• e*(p 3 )e{p 4 ). 
The hadronic cross section is given by 

da pp ^e*e{s) = ^2 J dr —^-da qij ^e*e{s), (75) 

q=u,d,c,s T ° 

where To = 4M| / s and the differential partonic cross section is 

2 dPS (2) 

d(t q q->e*e(s) = I M q q^e*e {§) | ■ (76) 

Here dPS*- 2 -* = dt/87rs is the two particle phase-space element and s = ts, where s is the hadronic center- 
of-mass energy squared. As it is well-known the parton luminosities are given by 



dr 



1 ~ J {fa/AiXitifb/B^Vi) + fa/B^, fj)fb/A(x, l4) > 



(77) 
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where the functions f a /A{x, n) are the particle distribution functions (PDFs). The amplitude squared for 
these processes can be written as 

\M q ^e(s)\ 2 = \ {Ut - Mf) \A q ^e*e{s)\\ (78) 

with s = (pi + P2) 2 , i = {pi - p 3 ) 2 , u = (pi- p A ) 2 and 

a ($\ — Cgg-yC^e ^C qq zCze*e C qq z BL C z BL e*e nQ , 

~ s + i-M*+ iMzTz + s - M\ Bh + iMz BL T ZBL ' ( 

where 

C qq -y = e q e, C 7 g1 \e L = e l e > Cze* T e L = ~ ^7, L e , (80) 

L L sin 2t)\y 

f f 

" smW Cw -sin2%' %^ " 6 /*/^ " ^L-y- (81) 

Here Lj = l| — e/ sin 2 and = — e/ sin 2 where l| is the isospin of the fermion f. Now, using 
the equations 

u = 2M?-i-s, (82) 



i = Ml-^ + yJs( S --Ml). (83s 



we can compute the cross section 

1 '■' dCT 



/i ri dL - 

dy dr-^ a{Me,y,T,s). (84) 



The numerical results for the selectron production cross sections are shown in Fig. 14 for different 



scenarios, with gsi assumed to be the maximum value allowed by the experimental constraints: Eq (28). 



We have compared our analytical results for the cross section with the results in Ref. [35 ] and found the 



same result in the case of the MSSM. In Fig. 14 we can see that even in the MSSM the cross section can be 
large and when the Zbl is included the cross section can be even larger due to the resonance enhancement. 
For example when the Mz BL = 1 TeV one can have a cross section above 1 fb where the selectron mass is 
below 450 GeV and y/s = 7 TeV. 

B. Signals with Multi-Leptons 

In this paper we wish to investigate the most promising signals associated with lepton number violation, 
through the process 

qq -> J,Z*,Z* BL -> e*e, -> efer^Xi -> efefefe* 4j, (85) 
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MSSM 

mz BL = 1 TeV 

- mz BL =1.5TeV 
m z =2TeV 



300 400 

WgJGeV] 
(a) 

FIG. 14: Drell-Yan production cross sections for the charged sleptons in our model. The dashed line corresponds to 
the prediction in the MSSM and the solid lines show the results in our model for different values of the Zql mass. 



The gauge coupling gsL is assumed to be at the maximum value allowed by the experimental constraints: Eq (28 1 



where i,j, k = 1..3 are generational indices. See Fig. [5]for the illustration of these signals. However, we 
will focus on the electron and muon channels since hadronic activity associated with the tau would blemish 
the signal of lepton number violation, namely three same-sign leptons, one lepton with opposite sign, four 



jets and no missing energy. Taking a cue from Eq. (38 1, which shows that the left-handed (right-handed) 
sleptons receive a negative (positive) contribution to their mass from the B — L D-term, we assume that 
only the left-handed sleptons are producible through this process. 

We begin by giving an estimate for the number of events in the limit of mostly bino, wino and Higgsino 
LSP. We present results for a 7 TeV LHC with 10 fb _1 of data. The combinatorics factor for the channels 
of interest are given by 

(86) 



F jk =2 (2 - 6 jk ) Br (ef -> efxi) x Br (xi -> efW^ 



x Br (xi -> 4 WT ) x Br ( w 



± 



so that the final states are 



e^e^e^e^Aj, e =c e + e + /x + 4j, e =c e + /x + ^ + 4j, 



J3, 



± C =F„T„=F, 



/i =t ^ =F e =F e =F 4j, / u =t /i =F e =F ^ =F 4j, /i ± ^ =F /i =F ^ =F 4j. 
We assume that Br (ef -»• efxf) ~ 100%, 33%, 100% for a bino, wino and Higgsino LSP respectively, 
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Hierarchy 


LSP 


Br(x? -> e^W*) 


Br(x? -> (t^W*) 


J~ ee 


**i eft 




NH 


Bino 


1-20% 


10-50% 


0.0001-0.036 


0.002-0.18 


0.01-0.22 


NH 


Wino 


1-20 % 


10-50% 


0.00001-0.004 


0.0002-0.02 


0.001-0.024 


NH 


Higgsino 


2-25 % 


10-50% 


0.0004-0.056 


0.004-0.22 


0.008-0.24 


IH 


Bino 


10-60% 


10-30 % 


0.009-0.32 


0.018-0.32 


0.009-0.08 


IH 


Wino 


10-60% 


13-30% 


0.001-0.035 


0.003-0.034 


0.002-0.009 


IH 


Higgsino 


10-60% 


13-35% 


0.008-0.32 


0.024-0.37 


0.016-0.11 



TABLE III: Ranges for the branching ratios of the LSP to charge lepton and W boson taken from the corresponding 
dense regions in Figs. 
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13 These are used to calculate the overall combinatorics factor, Fjk for the final state 



e i e T e j e k^3' Values are separated by the composition of the LSP: mostly bino, wino and Higgsino and for both the 
normal and inverted hierarchies. 



following the discussion in Section VII B The branching ratio of the W boson into jets is about 67%. For 



the RPV neutralino decays we pick the most prominent regions from Figs. 11-13 and display these values 
along with Fj^ in Table III Values are shown for both the normal and inverted hierarchies. 



To calculate the number of events expected after 10 fb 1 of data, we convolute the combinatorics in 



Table nu\ with the cross sections for a 1 TeV Zbl shown in Fig. 14 and multiply by ten; the results are dis 



played in Figs. 15-17 Notice that in most of the scenarios one can have several events which are basically 



background free. The main backgrounds coming from ttWZ and jjjjW^W^Z are very suppressed. 
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FIG. 15: Number of e^e?e^ej^4j at a 7 TeV LHC for 10 fb 1 of data for a bino LSP. Branching ratio values are 



shown in Table III while cross section values are taken from Fig 14 Data is divided into (a) for the NH, and (b) for 
the IN. 
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FIG. 16: Number of ef ef ef ef 4j at a 7 TeV LHC for 10 fbT 1 of data for a wino LSP. Branching ratio values are 



shown in Table III while cross section values are taken from Fig 14 Data is divided into (a) for the NH, and (b) for 
the IN. 
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FIG. 17: Number of efefefe^4j at a 7 TeV LHC for 10 fb _1 of data for a Hig gsino LSP. Branching ratio values are 



shown in Table III while cross section values are taken from Fig 14 Data is divided into (a) for the NH, and (b) for 
the IN. 



In order to understand the testability of the model at the LHC we show curves of constant number of 



e^efe^e^ij events per 10 fb 1 of data in Fig. 18 in the Br(x? —> e^W T )—Br(ef —> efxl) plane. Values 



200 GeV. In the case of the observation of 



are shown for a seven TeV LHC, with a 1 TeV Zbl and m e -, 
such events a the LHC, the Zbl mass can be reconstructed from electron-electron and muon-muon events 
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and the selectron mass may be reconstructible from its decay into two leptons and two jets. Therefore 

can 



allowing the calculation of the cross section for charged slepton pair-production. A plot such as Fig 1 8 



be used to get a better handle on the two unknown branching ratios and shed further light on the model. 
In order to estimate the reach of the LHC we also present curves of constant number of efefe^e^ij 



events per 10 fb _1 of data in the Br (ef — > ef Xi) ~ m e 4 plane in Fig. 19 This is again for a seven TeV 
LHC, with a 1 TeV Zbl and we show two possible values for Br ^x? ~~ * e fW^\ representing the upper 
(lower) part of that range in blue (green). One can see that even if the slepton mass is around 450 GeV one 
could observe a few events with multileptons and four jets. It is important to mention we satisfy the recent 
bounds coming from ATLAS Il36ll37l. 
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FIG. 18: Curves of constant number of events for the final state ef&fejej&j in the Br (x5 — > ejW^) 
efxi) plane. Values are shown for a seven TeV LHC, Mz BL = 1 TeV and = 200 GeV. 



Br (if^^- 



IX. SUMMARY AND OUTLOOK 

In this article we have studied in detail the theory proposed in Ref. 1 15 ] which we consider the simplest 
gauge theory for R-parity violation. This theory makes a prediction for the LHC since in order to break 
the B-L gauge symmetry the right-handed sneutrino must get a vacuum expectation value and one should 
observe lepton number violation at colliders. We have found the following results: 
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FIG. 19: Curves of constant number of events for the final state efefe^e^j in the Br (ef — > efxi) ~ m ei 
plane. Values are shown for a seven TeV LHC, M Zbl = 1 TeV and for two different values of Br (%i ~> ef\Y T ), 
representing the upper (lower) part of that range in blue (green). 

• In Fig. 2 we have illustrated in a simple way that one can have a realistic scenario for all sfermion 
masses even if some of the sfermion masses have a negative and large contribution from the B-L 
D-term in the theory. We have shown that in order to avoid tachyonic masses one should satisfy the 
condition Mz BL < \f2mj. This is a simple result which helps us to understand the constraints on 
the spectrum. 

• The full spectrum of the theory and the constraints coming from neutrino masses were analyzed in 
detail. The spectrum for neutrinos is interesting since it contains five light neutrinos: three active 
neutrinos and two sterile neutrinos. Using the experimental constraints on the masses and mixing 
for the active neutrinos we show in Fig. 4 the allowed values for the vacuum expectations of the 
left-handed neutrinos and the Yukawa couplings. As we have discussed in the text, these results are 
crucial to understand the decays of the lightest supersymmetric particle in the theory. 

• In Figs. 5 and 6 we have shown the properties of the new neutral gauge boson in the theory, the B-L 
gauge boson. Since one has two extra light neutrinos in the theory the invisible decay width is larger 
in this case. The contributions of the supersymmetric particle to the decay width are small and so the 
Zbl is like the B-L gauge boson in the non-SUSY scenarios. 

• We have investigated the neutralinos decays in great detail. In Figs. 7-9 we have shown the results 
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for the decay length in the different cases. As one can appreciate in Figs. 7-9 there are some sce- 
narios in the Bino limit where one could expects displaced vertices. The branching ratios have been 
investigated in Figs. 10-12 and we can summarize the results in the following way 

Br(x? -> rW), Br(x? -> fiW) > Br(x? -> eW), 
in the Normal Hierarchy and 

Br(x? -> eW), Br(x? -)• /iW) > Br(x? -> rW), 
in the Inverted Hierarchy, in majority of the parameter space. 

We have studied the main production channels for the charged sleptons at the LHC. In this case one 
can produce the charged sleptons through the photon and the Z as in the MSSM, and through the 
new neutral gauge boson, Zbl, in our model. As we have shown in Fig. 13 the production cross 
section can be large and thanks to the presence of the Zbl one can have even larger values for the 
cross section due to the resonance enhancement. We should point out that this production channel 
(throught the photon and Z) is very important to understand the signals in any model for R-parity 
violation. 

The most striking signals for lepton number violation in this context are the channels with three 



leptons with the same electric charge and four jets. In Fig. 15 19 we have shown that one can have 
a large number of events at the LHC with only 10 fb _1 . The background for these channels is 
suppressed, therefore there is a hope to test or rule out this theory in the near future. 
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Appendix A: Mass Matrices 



In the case of the CP-odd neutral scalars, in the basis (Al, Ar, Ad, A u ), one finds that the mass matrix 
reads as 
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B„ 



M 2 



odd 



m B, 



V2 



^B p 
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V2 
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V2 



a u v R 
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.J-Y nvz, Ly ,, VT «&R/, I YvV-vlvr 



V2 



a v v R 
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Bfi 
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(Al) 



while for the CP-even scalars, in the basis (hi, h R , h d , h u ), one finds: 



where 
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a uH — 



\ (a! + 9 2 2 + 9 2 BL ) v\ + %B V -\ (g% L - 2Y?) v L v R - B v \ 
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V -l(5i 2 + ff 2 2 )^-^ \{9l+9l)vl + ^B^-^j 
In the case of the charged scalars, in the basis (e, (e c )* , H~ , H^) the mass matrix reads as: 



Cfi B e 

Be Cf 2 



r 2 



(\glv L v d - \Y^v L v d - -^Y u fiv R \g\v L v u - \Y^v L v u - -^a v v R \ 
\Y e Y u v R v u + -^a e v L \Y e Y u v R v d + ^Y e /j,v L J 



5/i + \g\v u v d \g\ {y\ + 4) + M Bfl - \Y?v\ 



In the above equations C\ x and C 2 2 are given by 

/-(2 1 2 z' 2 2 \ i t^2 2 1 t^2 2 , U -R d 

On = 4S2 K - «d) + - 2^ v « + — B ^ 

and 

^ 2 2 = mi + ^ 2 (. 2 - «3 - ti) + (4 - 4) + \y 2 (vj + 4) . 

We also define for convenience: B v = 4= (Y u fiv d — a y v u ) and £? e = 4= (Y e /j,v u — a e v d ). 
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The amplitude for Zbl are 



Appendix B: Decay Amplitude 



2 4 / 9BL f \ 2 2 I ^ f \ 

M{Z B l ->■ fifi)\ = z c f\~^ n BL) m z BL I 1 + m 2 1 > ft = u,d,c,s,b,t,e,fj,,r; (Bl) 



\M(Z B L -> Witt = o ^T%L m Z Rr , 



\M(Z BL -»• ATiV)| 



M{Z BL ^f a fV 



3 V 2 

V^/k^J A _ 4 m ^ 
3 I 2 ™^ I 1 4 m| M 



1 / $BL / \ 2 2 



/ 2m 2 ; +2m 2 ? f 



2 2 

m - — 

fa fp 



\ 



Ly 7X1% 

^BL ^BL 
faf/3 = QiaQip, kal*p, , URi^Ri- (B5) 



(B2) 
(B3) 

(B4) 



Here, i is a generation index, c/ are color factors (c 9i = 3, = 1) and are the unitary sfermion mixing 
matrices. 

Appendix C: Neutrino-Neutralino Mixing Matrix 



In the basis ifj T = (u, x), the mass matrix has the general form 



M 



where Ai is diagonalized by J\f 



M ] MM* 



03x3 rn D 
ml M X/ 





M® , 



(CI) 



(C2) 



is the diagonal mass matrix for the light neutrinos, M x is the diagonal mass matrix for the neutralinos 



and 



TV 



U V 

v c u c , 



(C3) 



Eq. ([C2l» yields 



m D = U ] m D V* + V}m T D U* + V^M X V C , 
M° = V^m T D V: + U\M X U* C + Utm T D V\ 
= U ] m D U* + V}M X U* C + V^m^F*, 



(C4) 
(C5) 
(C6) 
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and 



E/~ U c ~ 0(1); 



,m 



D 



V ~ V c ~ 

m D ' 



(C7) 



The unitarity condition yields the following expressions 

uu ] + yy f = u ] u + y c f y c = v c v} + i/ c e4 = + u\u c = i, 
uv} + vut = u ] v + v}u c = 0. 



(C8) 



In Eq.|C6|the V^m^V* is negligible whereas in Eq. 



C5 



the M x term dominates. Therefore 



U ] m D + V}M X = 0, 

UlM x U c = M°, 
m D U* = VM, 



(C9) 
(CIO) 
(Cll) 



where the last expression is a result of inverting Eq. |C10[ substituting it into Eq. C9 and making use of the 
unitarity conditions. Substituted Eq. C9 into Eq. |C4| yields 



V}m T D = m°U T 



These results can be use to manipulate the seesaw relation: 



m v = Um^U T = rnoM 1 mJ, 



(C12) 



(C13) 



where m u is the nondiagonal light neutrino mass matrix diagonalized by U. Substituting Eq. (C 12 1 for 
m®U T , rearranging using the unitarity condition yields and solving for V yields 



V = m D M- x XJ c 



(C14) 



This can be rewritten by inverting Eq. (C5 1: 



V = m D U* (M. 



(C15) 



Where V can be identified with Vi a , the matrix that describes the mixing between the neutrinos and the 
neutralinos and is necessary for computing the neutralino decay properties. This result agrees with the 



naive expectation from the mass insertion approximation. While in the decay widths, Eq. (70 74), factors 
of U and E (the matrix that diagonalizes the charged lepton mass matrix) appear, they do so as sums of 
unitarity quantities and therefore are either zeroes or ones. 



35 



Appendix D: LSP Candidates and Their Final States 

The violation of R-parity increases the space of possible LSPs, which is now no longer restricted to 
chargeless fields. We therefore take the time here to make a quick survey of the possible final states in 
this model. For each possible LSP, we consider its production, if lepton number violation is observable 
in principle and if there are any obstruction to this observation. Finally, we judge which LSP leads to the 
most interesting signals. For us, these are the signals where lepton number violation is explicit: same-sign 
leptons with no missing energy (which might be due to neutrinos thereby confounding the counting of 
lepton number). Of course this can only arise from neutral LSPs. 

• Gluino (g) LSP: 

Gluino pairs are produced through strong cross sections at the LHC. Their possible decays are 

PP ~^ 99 ~ ttbbe^ej, tttivu, 

PP -> 99 -> 4j ef ef, 4j uu, 

where the former is favored if the third generation squarks are lighter than the first two. The gluino 
decay width can be estimated as 

i x+ 

For Mg = 100 GeV, Mq = 500 GeV, M%+ = 500 GeV and v\ = 10 MeV, one finds that the decay 
width is smaller than 10~ 13 GeV: a long enough lifetime for the gluino to form bound states but short 
enough so that it decays within the detector. 

In principle, these channels can yield spectacular signals at the LHC. However a recent inclusive 
analysis in the search for isolated same-sign muons published by the ATLAS collaboration has placed 
a model independent upper bound on the gluino pair production cross section of 58 fb [34]. Imposing 
this bound translate into a lower bound on the gluino mass of around 1 TeV indicating a heavy SUSY 
spectrum. Since we are interested in the scenarios with low energy supersymmetry we do not pursue 
this scenario further. 

• Squark (q) LSP: 

A stop LSP allows for final states with two third generation quarks of the same type and two leptons: 

pp — > t*t — > bb ef ej , or tt vv, 
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while a first or second generation squark LSP has channels with two jets and two leptons: 

pp ->• q*q ->• 2j ef ej , 2j i/i/. 

These channels have strong cross sections but do not provide information on the violation of the total 
lepton number. Since squarks in this case act as leptoquarks (each decaying into a quark and lepton) 
bounds on this scenario can be derived from leptoquark searches. 

• Charged slepton (e"j) LSP: 

Charged sleptons can be pair produced through the Z and Zbl, with signals 

pp — > e*ei — > it bb, efe^ w, 

where once more, lepton number violation is not discernible. The tb finals state is due to the mixing 
of the charged sleptons with the charged Higgs boson which typically decays in this way. The 
leptonic channel arises due to the R-parity violating mixing between the charged leptons (neutrinos) 
and the charginos (neutralinos). 

• Sneutrino (z/j) LSP: 

Sneutrino pair production also proceeds through the Z and Zb l with the following possible final 
states: 

pp — > v*v — > bbbb, ww, efejefe^ 

The first final state results from the R-parity violating mixing of the sneutrino with the Higgs boson, 
while the latter two states are due to the R-parity violating mixing between the charged leptons 
(neutrinos) and the charginos (neutralinos). 

• Chargino (x ± ) LSP: 

Charginos pair production is possible through the Z and leads to channels with two charged lep- 
tons due to the R-parity violating mixing between the charged leptons (neutrinos) and the charginos 
(neutralinos). 

pp ->■ x + X~ e t e j zz i wW + W~ . 
While in this case lepton flavor violation is observable, total lepton number cannot be probed. 
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Neutralino (x§) LSP: 

This scenario allows for several interesting channels with lepton number violation. If the neutralino 



is Higgsino-like the pair production (See Section VII C for information on neutralino decays): 



PP -> X1X1 -> 4 J ef e^, 

is possible through the Z, as well as associated production which gives rise to channels with three 
charged leptons: 

PP -> Xixf -> 4? ^ ef . 

Unfortunately, these channels are interesting only in the Higgsino-like neutralino scenario and in 
general the cross sections can be small. However, striking channels with three same-sign charged 
leptons, multijets and no missing energy through the pair production of selectrons are generally 
present: 

PP -> e*ii -)■ ef ef ej ef 4j. 

Such striking signals maybe the signatures that help test this model at the LHC. As we show above, 
the production cross section can be large and there are no relevant backgrounds. 
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